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\0 ■ Abstract 
O ; 

We demonstrate how to resum phase space logarithms in the Sterman- Weinberg (SW) dijet decay 

(N . 

' rate within the context of Soft Collinear Effective theory (SCET). An operator basis corresponding 

to two and three jet events is defined in SCET and renormalized. We obtain the RGE of the two 
and three jet operators and run the operators from the scale u 2 = Q 2 to the phase space scale 

| u 2 = 5 2 Q 2 . This phase space scale, where 5 is the cone half angle of the jet, defines the angular 

region of the jet. At (J, 2 we determine the mixing of the three and two jet operators. We combine 

. these results with the running of the two jet shape function, which we run down to an energy cut 

' scale fig. This defines the resumed SW dijet decay rate in the context of SCET. The approach 

\o . 

outlined here demonstrates how to establish a jet definition in the context of SCET. This allows a 

program of systematically improving the theoretical precision of jet phenomenology to be carried 

Qh! out. 
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I. INTRODUCTION 



The phenomenology of QCD lets has been examined over the years with increasingly 
sophisticated techniques , 0]. Recently, an effective field theory (EFT) of QCD 

containing collinear degrees of freedom interacting with soft and ultrasoft gluons, SCET; jj], 



was formulated. Using SCET, nonperturbative corrections to dijet event shapes have been 



examined in a series of papers 
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In this paper, we use SCET to sum large phase space logarithms 1 present in the pertur- 
bative expansion of the dijet decay of the Z boson. Our approach is general and addresses 
the problem of establishing a jet definition within SCET. We have chosen our initial state 
to be a Z boson to avoid complications from QCD interactions with the initial state. The 
logarithms resummed have the form a" log n+1 in fixed order perturbation theory. At 0(a s ) 
we are resumming the a s log (2/3) log (8) double log in terms of the jet parameters 8 and 
(3. We also resum a class of a™ log™ logarithms given by a™ log™ (§) terms in the fixed 
order perturbative expansion. Here 8 and (3 are parameters that impose phase space cuts 
according to the SW jet definition. 2 

The SW jet definition defines the dijet decay rate as an integration of the triple differential 
decay rate over a phase space restricted by cuts. The cuts on the energy and angular 
integrations define the dijet decay rate to be the sum of three partonic final states. The first 
state (SW1), is comprised of a quark and anti-quark contained within different jets, where 
the jets are angularly defined as cones of half angle 8. The second state (SW2), is comprised 
of a quark and anti-quark each contained within a different jet, plus a gluon (unrestricted 
in direction) with energy E g < (3M Z . The final state (SW3), is defined as a quark and an 
anti-quark within different jets and a gluon of energy E g > (3 M z within the jet cones. 3 

In this paper, we match the amplitude for the decay of the Z onto an operator basis 
defined in SCET. This defines contributions of two and three jet final states at the scale 
Mj. We run down from the scale M% to the scale defining the dijet through the angular 



1 It has been shown that phase space logs can to be resumed for the decay B — > X c £is in 12 1. 

2 The SW jet definition is not the only possible definition of a jet. In fact, SW jets are somewhat experi- 
mentally disfavored compared to the k± algorithm 0, Q| ■ The relative theoretical simplicity of the SW 
jet definition motivates us to demonstrate the origin of the relevant logarithms using this definition. 

3 Note that both 8 and (3 are small (<S 1) in this definition; however one must also have (3 less than 5 such 



that sm{6) > /3/(l - /3)[l5] 
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cut on phase space. This cut scale is given by u| = S 2 M|. The usoft degrees of freedom 



are matched onto a dijet shape function as in |8|, 10, 11]. This dijet shape function is 
restricted to emit gluons that take away energy E g < (3M Z . This requires the introduction 
of a further phase space scale corresponding to this energy cut. This cut scale is given by 
[Ap = <5 4 Mf /£> 1_<5 , where B is an 0{1) parameter. We fix B by matching onto the SW dijet 
decay rate after we run down to the jet scales. We discuss the physical significance of these 
jet scales in Section III] in more detail. 

?he resummation of Sudakov [ijj] logs for jets has been investigated in Q, Q, Q, 
2l| for various jet definitions. The leading logarithm results reported here have been 
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produced in these other contexts. We approach the re-summation with the EFT formalism 
as this approach offers several advantages. For example, the EFT construction includes 
nonperturbative effects through the matrix elements of the usoft degrees of freedom. Within 
the effective theory, one can calculate to high orders in RGE improved perturbation theory 
in a clearly defined and easy manner. As an example of the utility of the effective theory 
approach, we have determined the mixing of the three and two jet operators. This mixing 
determines a class of NLL contributions in the RGE improved perturbative expansion. 

It is important to note that although the SCET final states discussed here are an inclusive 
sum over the set of hadronic final states, the sum is not a colour singlet. Gluons are 
exchanged between the low energy partons comprising the jets in the hadronization process. 
This is a difference in calculating dijet production compared to the use of the operator 
product expansion in other inclusive decays, such as semileptonic B — > X c iv decay 22]. 
Experimental results 23( indicate that this issue does not invalidate the SW jet definition. 

The outline of this paper is as follows. In the remaining parts of the introduction, we 
calculate the 0(a s ) phase space logarithms of the SW jet definition. We then briefly review 
SCET. In Section [XT] we formulate the phase space cut scales corresponding to the SW jet 
definition. In Section ITO1 we discuss the matching of T(Z — » qq) onto the dijet operator. We 
renormalize the operator basis and run down to the the scale /j, 2 -. We match the contribution 
of collinear gluon emission in Z — »• qq g onto three jet operators in Section llVl and renormalize 
the three jet operator basis. Section IIVBI is devoted to the running of these operators. In 
Section |V] the mixing of the three and two jet operators is discussed. In Section IVHI we 
introduce the 'Phase Space Renormalization Group' and run to the phase space cut scales. 
Finally, we obtain the leading log resumed dijet decay rate in SCET. 
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A. Phase space Logarithms at 0{a s ) 



First, we review how the phase space logarithms present in the QCD calculation of 
Y[Z — > J J) are generated by the SW jet definition. Recall that to 0(a s ) the amplitude 
for F(Z —> J J) is given by the following expression, grouped in terms is qq and qqg final 

states 



A 2 



\ A + Av + A W i + A W 2 \ + I Am + A B 2 \ R ■ (1) 

The subscript R refers to the restricted set of qqg final states in the jet definition. The 
amplitudes above are the tree level amplitude A®, the vertex correction amplitude Ay and 
the wavefunction amplitudes Awi- The Am are the bremstraulung amplitudes. 

Consider first Aq and the corrections Av, Awi, Awi- These contributions to T(Z — > qq) 
can be determined by calculating the imaginary part of the forward scattering diagrams in 
Fig. (1). We calculate in the massless quark limit throughout this paper in the rest frame 



Im 




FIG. 1: The imaginary part of the forward scattering amplitude in QCD gives the qq final state. 
The operator insertion is the weak neutral current j nc = qT a qe a with T"" = 7 " (gy + g& 75). 
The wavefunction renormalizaton graphs are not drawn as Landau gauge can be used so that the 
wavefunction contributions vanish. The sum of imaginary part of the diagrams is gauge independent 
at 0(a s ). 



of the Z , and always in d = 4 — 2 e dimensions. We average over the initial polarization of 
the Z in d — 1 dimensions, and obtain the decay rate 



r(z-> g g)(/0 



ff2 -(4 + ^)[Mr £ (4,) 2 ^fi 3 - 2£ ]C(/,). 



(2) 



We have presented our results in the form suggested in {9] and agree with their result. The 
coefficient in MS is 



cfi^) = i + as ^ )CF 



7T 



|)(.n[^M]_I ) _ 4 + 

— LL c 



1 i„2r- 2 2VP?l 
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In 



0(a*)(3) 



The qqg final state amplitude is Az mg —\ Abi + Ab2 | 2 and is obtained by calculating the 
bremstraulung graphs. The double differential rate is 




FIG. 2: The three body qqg contributions. 



dT 1?9 _ N C C F {g 2 y + g|) 1-2 £ (4 7T) 4e fl 3 _ 2 e Vt 2 _ 2 e (d _ 4) (d-4) ( _ 2) (e) pr , ^ 

cte 9 dx q ~ 256 vr4 (3 - 2 e) J ^ ^ Z " ^ j 

where the energies of the decay product quarks are x q = 2E q /M z and x q = 2E q /M z . The 
angle between the quarks has been fixed by a delta function to be z = 1+2(1— x q —x q ) / (x q x q ). 
The function F is kept to order e 2 and is given by 



Integrating this expression over the remaining three body phase space in d dimensions gives 
r(z - = + ^)[M|- 2e (47r) 2e |^n 3 _ 2e ]^(/,), (6) 

with 



ego,, = ((I + 2) (lD[ M, _ I, _ » + £ _ i ^f) + 0(aS . (7) 

The full decay rate to is the result 

12 7T y 4 7T y 

The SW jet definition introduces cuts on the energy and angular integrations of Eq 
(jl]). The resulting phase space integrations are difficult, some details can be found in {{J. 
Retaining the terms that diverge as the energy fraction f3 and the cone angle S go to zero, 
one obtains 



C™"»(M) = ((I + |)(l n [Ml] --)-- In 2 [M] ] 



(_ 4 ln(2 3) ,„(,) _ 3 , nW + f - i£) . (10) 



Combining this result with Eq.([51h one obtains the SW result; we find 

F(Z 7J)(M») = N cMJ& + &) f 1 + «^ 8 _ £ _ 3In4 _ \ 

VZtx \ n I 6 J 

Note that the result is free of explicit IR singularities. This is expected as the SW jet 



definition is constructed to satisfy the KLN theorem 32j, [33[- However, the result does 
depend on logs of the cuts used to partition phase space. 

We will now examine this decay in SCET. The different final state phase space configu- 
rations will be separated out into different operators using the invariant p g -pi- By matching 
and running, the renormalization group allows us to sum the Sudakov logarithms and define 
a dijet amplitude. 

B. Jets in SCET 

SCET contains collinear and ultrasoft degrees of freedom that are relevant for QCD jet 
studies. 4 Consider a final state gluon produced by the decay products of the Z. The gluon is 
collinear in direction when it has scaling p^ = Mz(X 2 , 1, A), where A is a small parameter 
in the effective theory. We have used the lightcone coordinates L M = L~ , L±) such that 
L + = rti ■ L,L~ — T%i ■ L and = L 11 — L + nf/2 — L~ /2. The lightcone vectors used 
here are defined as = (l,nj) and = (1, — Hi) with rij a unit vector. 

Alternatively, the gluon could have ultrasoft scaling p us = M Z (X 2 , A 2 , A 2 ). When A ~ 



Aqqd/Mz, one cannot perturbatively expand the decay rate in terms of usoft gluons. The 
nonperturbative effects of these gluons are given by the matrix elements of ultrasoft degrees 
of freedom, which give nonperturbative corrections to the decay rate 8|, |9|, [l0|, For this 
reason, we take our final state gluon radiation to be collinear in our perturbative expansion. 
Expanding the decay current j M in terms of collinear and ultrasoft degrees of freedom, the 
SCET jet fields are given by 

Xn q = [£,n q W nq \, 

Xn q = [W^n 9 ]- (12) 



See l7f for more detailed reviews. 
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The £ field is defined with the large momenta p in the direction n« removed 



using the projector 



V 



P„ 



Also present to preserve collinear gauge invariance are Wilson lines, 



W ni (q) 



E ex P [~m ni ' An ^i 

.perms ^ ' 



(13) 



(14) 



(15) 



E ex P ■ y 



.perms 



(16) 



II. THE SCALES OF THE SW CUTS 



Three small parameters exist in the construction of SW jets Aqcd /Mz, 8 and (3. We take 



5 ^> J Aqcd /Mz for the cuts to be above the hadronization scale. To sum large phase space 
logs with the renormalization group, we must associate the phase space cuts with scales n 2 cut . 

For the angular region defining the jet (5), we utilize the approach of 2], [], [24J] and take 
a collinear momenta comprising the jet to have collinear scaling pj" 1 -* ~ Mz{5 2 , 1,5). This 
gives the cut scale of 



l4 = FMl. 



(17) 



We set the invariant mass of the initial off shell quark and anti-quark to this scale, P 2 2 — 
5 2 M§. Physically this scaling has a simple interpretation. The jet cone is defined by the 
half angle 5. For all momenta making up the jet,when z is taken as the axis of the jet cone, 

|P±I 



tan (S) 



(18) 



The usoft gluons radiate into the final state and deposit a small fraction of center of mass 
energy, (3 Mz, isotropically in the detector. For these gluons, we utilize the fact that this 
radiation cannot disturb the back to back configuration of the jet cones. The sum of usoft 
gluon emissions is restricted to carry a total energy (P°) out of the jet cones that satisfies 
sin(5)>P u 7(l-P°) 



151 ] . This requires the momenta emitted outside the jet cones to have 



scaling P u ~ M z (8 2 , 5 2 , 5 2 ). We take the energy cut scale of the usoft degrees of freedom to 
be 

^ = Mi' (19) 
where B is an 0{\) parameter. The correct resummation of leading logarithms in the SW 
jet definition requires that nfi/f^ 2 oc 5 2 . Physically, this corresponds to the usoft radiation 
not being able to cause the quark and anti-quark to lie outside the jet cones. Note that the 
results of the SW phase space integrals in Eq. ffTTj) are expanded in S and (3. All dependence 
on the jet parameters 5 and (3 that does not diverge as j3, 5 — > is neglected in Eq. ffTTj) . In 
a similar manner, the dependence on £>, 5 in our RGE is expanded. Dependence on these 
parameters that does not lead to large logs as j3, 5 — > is neglected. In this manner, any 
/x| dependence on B that gives the same retained contribution to the RGE will lead to the 
same double Sudakov logarithm. The quark and anti-quark B dependence is sub leading 
in the parameter 5 and neglected except for the component that scales as 5 2 . How the 
particular components of the usoft gluons depend on B is fixed when we match onto the SW 
jet definition. We will demonstrate how these phase space scales will establish the SW jet 
definition in SCET in Section IVTl 



III. Z^qq IN SCET 



At the scale fi 2 = M 2 we match T(Z — > qq) in QCD onto the dijet operator 25| 

0° 2 = XnJ* Xn r (20) 

The EFT loops are scaleless and vanish in dimensional reqularization at 0(a s ). The Wilson 
coefficient C 2 and the remormalization factor Z 2 de pen d on the large label momenta com- 
ponents 2p q ■ pg ~ M| of the quark and anti-quark 23]. The Wilson coefficient is given by 
the one loop matching condition 

((%)(!*) C 2 Qm) 



determining 

CM 



i + 



<q\A\Z 9 ){ii) 



a s (fi 2 )C F/ 3 l -2p q -p q - 



(21) 



2tt 



-In f- 



7T 



1 | a a (fx 2 )C F 1 | 3 



2vr 



2r 



In 



4 + ln z 

12 2 

-2p Q • Pq 



-2Pq ■ Pq 



}) + 0(aj). 



(22) 



The contribution from this operator at the scale /i 2 = M§ is given by Eq. (I21~T) with 

C 2( M z ) = l + a -^L{,i- i+ 7 ^y 0( al). (23) 

Note that contribution to the decay rate from the imaginary part of the Wilson coefficient 
vanishes at 0(a s ). The imaginary term can contribute at 0(a 2 ). 

To determine the counter terms in the effective theory, we use offshellness to regulate 
the IR. Both soft and collinear loops contribute to the UV divergences that determine the 
counter term. If one adopts two stage running 5 , the subtraction points of the usoft and 
collinear loops are taken to be the same, ie fi 2 = fi 2 . With this choice, the anomalous 



dimension of O2 is given by 



25] 



72 = -2a s 



da s 

«■<"> C F (log + a, + a(al) , (24) 



Note that we have used the MS result [26j for the epsilon poles of the counter term Zf . 

The Wilson coefficient of the operator at this lower scale is obtained by solving the RGE, 
with (3q = 11 — 2nf/3. For two stage running we find 

a < m= [mmT)) [mI) c ^ h) - (25) 

IV. Z^ggc/INSCET 

Beyond tree level, we match the contribution of T(Z — * qqg) onto the two and three 
jet operator basis. The QCD amplitude is matched onto a two jet operator when the 
collinear gluon emission can be described by Cscet and P 2 9 < 5 2 Mf or P? g < 5 2 Mf , where 
= P" + Pg. Conversely, when Pf g > 5 2 Mf for both i — q,q, the collinear gluon emission 
is matched onto a three jet operator. At the scale /1 2 = M 2 , we have the following amplitude 
depicted in Fig. (2) 



Abi + A B 2 = u(p q 



1-2 



-^%n (iM r» - (iM (-^ 7 ,T Q 

V gg / V qg/ 



v(p q ) s %e v . (26) 



5 We find it convenient to adopt one stage running in Section IVIl However we utilize two stage running 
initially to demonstrate the benefit of one stage running in Section IVT1 
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For the three jet operators, we take 2p g ■ pi, ~ M| and the intermediate propagator is not 



(27) 



resolved between the scales Mf and 5 2 Mf . Our matching to order A°Mf is given by 

A B1 +A B2 = J2 ^( q qA^\Of\e a )^), 

i=3,4,5 A iW 

where we have suppressed the summations over the light cone directions n^. For different 
phase space configurations, 2p g ■ p^ takes on different values between Mf and S 2 Mf. In 
running between Mf and 5 2 Mf the logs generated by the separation of scales depend on 
the phase space configuration of the event. Our formalism takes this into account by rii 
dependence in the anomalous dimensions of the three jet operators. 

To decompose the amplitude into operators with field content defined in SCET, we first 
rearrange the three parton amplitude into the form 



A 



V+A B 2 = U S2 (p g ) 

+ u s *{p q ) 
+ u s \p q ) 



(P g ■ P q )(P g ■ P 9 ) 

Pg ' [Pg ~ Pq)(Pgg e % ~ e 1 Pgp 

2(P g -P q )(P g -P q ) 
P g .(P ri -P g ) ie ^^e a a e' a P g ^ 5 
2(P g .P q )(P g .P q ) 



r°(g a T a )v*(p 9 ) e' a 

Tp{g s T a )v^{p q ) e' a 
r v (g s T a )v^(p q ). 



(28) 



We now expand such that the partons are collinear particles in the directions rii, ti 2i ti 3 for 
p q , Pq and p g respectively. The operator basis in SCET is dictated by the constraint of 
collinear gauge invariance in the three directions n 1; n 2 , n 3 . We find the following operators, 



^3 



°4 



<> 



Xn q Cf(n 3 ■ p 3 , m) [Wl g s W n3 ] T° X n< 
Xn q C A {ui -Pi, [Wl 3 g s G^f W n3 ] T Q Xn 9 i 
Xn q C 5 (ni -pi, Hi) [Wl g s 75 W n3 ] T v x n 



(29) 



4 Tli n 2 



The Wilson coefficients are of the form 

C 3 /3 (^n 3 -p 3 ,n i ) - 

C^i^Ki - Pi, Hi) -- 



ni ■ n 3 ) (n 2 ■ n 3 )(n 3 ■ p 3 ) 2 
2 



(30) 



\ (m • n 3 ) (m • Pi) (n 3 ■ p 3 ) (n 2 • n 3 ) (n 2 ■ p 2 ) (n 3 • p 3 ) 
where the non label part of the Wilson coefficient is 



(31) 
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with Cf = 1. We have used the O(X ) collinear gluon field strength g s G%f = [iD™ + 
g s A^ 3 Pg ,iV^ 3 + g s -^n3,pg\ an d &nz = i ^ a 13(7 11 G /3 * s the corresponding dual field. The am- 
plitude expressed in terms of these operators requires a summation over lightcone directions. 
This summation is such that when rii and rij label distinct collinear directions, they lie in 
different jet cones. This constraint along with momentum conservation is 

Csi/J 2 ) = CsO 2 ) [rii ■ nj - (1 + cos 25)] 8{n x ■ n 2 + n x ■ n 3 + n 2 ■ n 3 - 2). (32) 

We suppress these summations and constraints until Section IVl 

The three jet operators contain the same field content. In SCET these operators are 
expected to have the same renormalization 30|, |34|]. We present the operators in the form 
of Eq. (j29|) as they have a clear field content and Dirac structure. We also find this form 
convenient for the calculations in Section |V] 

A collinear gluon with large angular separation from both the quark and anti-quark can 
also be produced when the quark and anti-quark are collinear. In the original calculation of 
SW, these states did not contribute at 0(a s ). In the effective field theory, the corresponding 
operators vanish at leading order in A. This is another example of the the utility of SCET. 



A. Renormalization of the 3-Jet Operators 



The three jet operators are renormalized by determining the UV divergences in SCET 
using the Feynman rules stated in the Appendix. The UV divergences are regulated using 
dimensional regularization and MS. We utilize offshellness to separate out the IR diver- 
gences. 6 For some related details on one gluon external state renormalization calculations, 
27]. We also note that we use the background field method 0, Sj- 



sec 



As we examine perpendicular polarized collinear gluons in the n 3 direction, this causes 
the contributions from a number of diagrams to vanish in the effective theory. The remaining 
diagrams are shown in Figure 3. 

Calculating in the rest frame of the Z, the perturbative expansion of the operator can be 
determined. We obtain the following for diagrams (a) and (b) 



8tt 2 



1 



log 



-P 2 

q 

» 2 



We have used Feynman gauge. 
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(c) 



(d) 





FIG. 3: One gluon external state diagrams for the operators Oi in the effective theory. Collinear 
gluons are drawn as springs with lines, usoft gluons as springs. 



g 2 s C F C™ 



-log 2 

2 & 



» 2 



+ log 



(UUA^\O°\0) b 



g 2 s C F C™ 



1 1 1 

-— - — + — log 



8tt 2 



--log 2 
2 6 



uv L uv 
2" 



-Pi 

J 2 



log 



-II 

A* 2 



-p? 

q 

VP 

-P 2 

g 

A* 2 



71 

2 + — 
12 



(33) 



2 + 



71" 

12 



(34) 



We have defined C va = (£ ni £ n2 A^ u \0% |0) in the above. Defining the following function 



KftP^mliPjw)] = log 



-2Pi-Pj 



log 



-P 2 

3 

n 2 



log 



-log 



1 nj ■ Wj / rij ■ PA / PJ 



\ni-pj \rf 



(35) 



the loop results for (c, d, e) can be expressed as functions of a master loop integral that gives 
the following form 



/ 1 [(P i) n i ) ) (P ij n i )] 



1 1 1 

2tv) k 2 + ie P 2 /{ni -P^-m-k + ie P 2 /(nj ■ P,) -n y k + ie 

12 



1 



7T^ Tii ■ Tl-i 



-K [(Pi,ni), (Pj,nj 



^K[(P i ,n i ),(P j ,n j )} 2 + ^ 



Diagrams (c) and (d) give the following results: 

g 2 C A C™n 2 -n 3 



(UUA^\O°\0) c 

(UUAif\o°\o) d 



g 2 s C A C™ ni -n 3 



h[(P g ,n 3 ),{P q ,n 2 )}, 
h[(P g: n 3 ),(P q ,n l )]. 



(36) 



In obtaining these results, one must properly account for the zero bin to avoid introducing IR 
poles into the anomalous dimension. For this purpose, one can utilize the method outlined 
in 



291 ] . Diagram (e) can be calculated in a manner similar to 



where the anomalous 



dimension of was determined. We obtain 



(UL 2 A^\O°\0) e = -g 2 s (C F -C A /2)C^ ni - n 2 Ii [(P,,^), (P„ni)]. (37) 

The remaining diagram (f) can be determined by utilizing the full background field three 
gluon vertex |28j] and the two gluon feynman rule for the operator. For this diagram one 
finds 



9 2 C A C^ 

16 7T 2 

9 2 a C A C v ° 

16 7T 2 



log 



r_p2- 

9 


+ log 2 


r_p2- 

a 


TV 2 ' 


I ^ . 


I ^ . 


6 _ 



-2 



log 



-P 2 9 



I' 1 



(38) 



The remaining graphs are wave function graphs. For a soft gluon on the external quark 
field, the wavefunction graph vanishes as nf = 0. The collinear wavefunction renormalization 

Q 

graphs [7|] must also be added to this result. For each effective field £, with momenta Pi, we 
have the wavefunction renormalization factor 



Ait 



log 



A* 2 



+ 1 



(39) 



Combining these results and adding the corresponding counter term leads to the renor- 



malized perturbative expansion to 0(a s ). The renormalized operator is given by O 

)(0) 



(r) 



Z q J ZqO\ (Z 3 . As we utilize the background field method, the combination g A g is not 



renormalized 

ZM = i + 

4 



281 ] . The renormalization factor of the operator is given by 
a s (n) C F 



Air 

a s (fi) C A 



"2 3 2 / 


r p2 - 




r p- 2 i 


- + log 


i 


+ log 




e z e e \ 






» 2 _ 



All 

a s (fj,) C A 



-K[(P q , ni ),(Pg,n 2 ) 



K[{P q , ni ),{P- qi n 2 )\ 



-K [(P q , n 2 ), (P g , n 3 )] - -K [(P„ m), (P g , n 3 )] 



Air 



log 



p2H 
9 

V 2 



(40) 
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We scale the results by the one gluon matrix element (0 3 } = (£ ni £ n2 An^\(O%y r A 0) /C utT 
The perturbative expansion of the operator in the effective theory is 



All 

a s (jj) C F 



7- 



5tt 2 



3, 

2 log 



r pr 


3, 


r p2- 




r p 2 ~ 




r p2- 




q 

t?_ 


"2 


q 

t?_ 


+ log 2 


<? 

/i 2 . 


+ log 2 


9 

t?_ 





+ 



47T 

a s (/i) Ca 



K[{P q , ni ),{P q ,n 2 )\ 



a s (n)C A ( 5 7r 2 



47T 



12 



log 



p2 

9 

» 2 



+ log 2 



p2 

g 



8 



(-K [(P q , m), (P 9 , n 3 )] 2 - A' [(P ff , n 2 ), (P s , n 3 )] 2 + A n 2 ), (P„ m)] 2 ) (41) 

7T x 7 



The anomalous dimension can be determined from the renormalization factor. The scaling 
for the inner product p • Pj in this operator basis was defined to be 2P • Pj ~ Mf . Using 
this scaling, and taking the quark and anti-quark lightcone vectors along the jet direction, 

[rii ■ P^ / (rij ■ Pj). Using these results, we determine 

- L[(P h ni), (P i ,n i )],(42) 



we have (rij ■ Pj) / fa ■ Pi) (Pf/P, 



K[(Pi,ni), (Pj,rij)} 
where 



log 







r p- 2 " 




r p 2 i 


- log 


q 


- log 


q 




^\ 





L[(Pi,7ii), (Pj,rij)] = log 



2 rtj ■ Pj + rii ■ rij (rii ■ Pi 
2 (rii ' rij) n j ' Pj 



(43) 



The anomalous dimension is determined for two stage running to be 



Haiti, rii 



a s{n) C F 
Air 

An 



4 log 



Ml 



+ 6 + 4L[(P 1 ,n 1 ),(P 2 ,n 2 )] 





2 log 




+ 2 


Air 


[Ml_ 





(L[(P 1; m), (P 3 , n 3 )] + L[(P 2 , n 2 ), (P 3 , n 3 )] - L[(P 1; m), (P 2 , n 2 )]) . (44) 



Calculating the renormalization of the operators O4, O5 in a similar manner, one obtains 
the same result without any mixing complications, ie 7 3 = 74 = 75- The contribution of each 
field to the anomalous dimension is the Casimir invariant of the field's representation times 
the factor — a K (p ) log |/x 2 /M|| / (2tt). This result for the renormalization scale dependence 
agrees with 2jJ, [30(j. The remaining logs depend on labels which give the phase space 
configuration of the operator. 



B. Running the Three Jet Operators 



We have matched onto the operators 3 , 4 , 5 at the scale /1 2 = M|. One can run the 
Wilson coefficients of the three jet operators using the standard result 



Cj(fi 2 ) 

cm 



cxp 



HI H 



(45) 



14 



We retain the logs dependent on rij when solving the RGE as they are 0(1)7 Consider the 
logarithms dependent on rij. Expressing the inner products in terms of the angle between 
the vectors one finds a divergent log as 9 12 — > 0, 



L[(P 1 ,n 1 ),(P 2 ,n 2 )] =log 



2n y P j 



[l + cos o i:j ) (m ■ Pi) 



(46) 



2(1— cos 6ij) rij ■ Pj 

This limit is not of concern for the C F phase space logs, as when the gluon becomes collinear 
with the quark or anti-quark the fermion jets are back to back. The Ca phase space logs do 
diverge in this limit. 

We find that the Wilson coefficients Cj, for two stage running, are 

Oj^/i, Hi 



00 



00 



a 8 (M. 



Ml 



Ml 



Ci(M z 



(47) 



where 
Vx(M z ,m 



C F (3 + 2L[(P 1 ,n 1 ),(P 2 ,n 2 )\) - 



8tt(C f + C a /2) 



/WMf) 

+C A (1 + L[(Px, n x ), (P 3 , n 3 )] + L[(P 2 , n 2 ), (P 3 , n 3 )] - L[(P 1} m), (P 2 , n 2 )]) 



(48) 



V. THE MIXING OF 2 AND 3 

As one runs the jet operators down to a lower scale, eventually the phase space scale 
defining the jet cone is reached. When one runs below this scale, a third collinear direction 
can be resolved in the dijet. Consequently the two and three jet operators mix beginning at 
the scale /j/g . The following diagrams are defined in the effective theory by taking the limit 
before the loop integrals are performed. For the interacting collinear directions (rij, rij) this 
limit is given by 

lim A (rij • nj,rii ■ rij, ■ rii/j) = lim A (1 + cos (26 + <p) , 1 — cos {25 + <p) , 1 — cos (ir — § — 0/2)) 

where <p 5. This determines the mixing of the two and three jet operators at leading 
order in A and fixes a SW jet cone in the effective theory. We find the following mixing for 
3 and 2 



7 We thank A. Manohar for conversations on this point. 
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(9) 



(h) 



FIG. 4: The mixing of the three and two jet operators. 
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(49) 



Diagram (g) is given by this result with P q — > Pq. Note that these results agree with the 
form obtained in Eq. (1331) . Implementing the zero bin subtractions in these integrals ensures 
that the divergences present are all UV 29j. The mixing is given by the sum of (g, h) 



(tnJf»\OZ\0M)C F g* (2 2 fl 
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• (50) 



frl 

This renormalized matrix element is defined by (£ m £ n2 IO3 10) being free of UV divergences 
where 



z 



3,2 



(UUIO^IO)^)^^ 2 /2 3 /l 3 



8tt 2 Z 



3,2 



e 2 2e Ve 4 



87T 2 Z 3 ,2 V 2 6 2 S 

The matrix element requires the following counter term 



3\ / 




- loe 


4/ V 




r-p 2 " 
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. t4 . 





-P 2 

g 



+ - lot 



log 



_p21 
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(51) 



^3,2(^5) 



2tt 



/ 2 3 1 / 


\-PT 












+ log 






b + 2^-er g 


[ . 


[ t4 . 





(52) 



Taking P 2 = Pjj = S 2 M§ we find the off diagonal term of the anomalous dimension matrix, 



73,2(^5) 



3a s (/j, s ) Cj 
2~n 



(53) 



1(3 



For O4 and O5, the mixing of these operators into the two jet operator vanishes at 0(a s ). 
This occurs due to the operator Dirac structure, demonstrating the utility of simplifying the 
Dirac structure of the three jet operators. 

We can now determine the renormalized contribution of the three parton amplitude to 
the dijet decay rate. The contribution is the finite real part of Eq. dSlj) . given by 

(u UiOsfWrf) = (u uom cm) asi fl GF (l-^f)- ^ 

VI. PHASE SPACE RENORMALIZATION GROUP 

To determine T(Z — > J J)(fi$) we integrate over the two body phase space in d dimensions. 
The tree level result is 

= / (O^ZfO T?\ fO^fO JP.\ ^ (( 2)) 2 



(2^)3(2^) (2 7r)3(2E 2 ) 



states, pol 



62 71 A — Ze 

= ^(9 2 v + 9 2 A) + 0(e), (55) 

IZ 71 

so that 

T(Z JJM> = + A) (CTW) + C 3(A ) (f (56) 

Reexpanding this expression into fixed order perturbation theory, neglecting non logarithmic 
0(a s (Mz)) terms, gives 

T(Z -* J J)(Ml) = + A) (l + ^AMML (_ 3 log(i) - 2 log 2 W )) . (57) 

Comparing this result to Eq. (ITT|) . we see that we produce the purely collinear contribution 
—3 log [5] correctly, however the double Sudakov log in not correctly reproduced. The reason 
for this is that only part of the double log is present. We cannot run the collinear degrees 
of freedom of the dijet operator below the /x| scale. However, the usoft degrees of freedom 
run below this scale and are restricted by the SW cuts to lie below the scale These 
scales are related by /i| oc 5 2 /i|. To obtain the correct fixed order perturbative expansion 
for T(Z — > J J), this restriction of the usoft degrees of freedom must be included. Including 
this restriction, we find the correct expression to be given by Eq. (l72l . 
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To determine the correct result, we implement these constraints with a one stage running 
formalism using a Phase Space Renormalization Group (PSRG). 8 Running in this one stage 
formalism, the collinear degrees of freedom of the dijet run to the scale /x 2 ; simultaneously, 
the usoft degrees of freedom run to the scale /j/^. 



The usoft matrix element for the dijet is defined in [8|, |24|. The matrix element depends 
on the momenta k of the usoft degrees of freedom where 

S{k) 



£K*«(*)|oj|o>l 5 



(58) 



We have indicated the required renormalization of the dijet shape function. The usoft 



radiation has momenta k = Mz(\ 2 , A 2 , A 2 ) where A ~ J Aqcd/Mz and the cuts are taken so 
that 5 ^> Aqcd/^z- The cuts being larger than a typical component of an usoft momentum, 
in the effective theory, no cut restriction 9 is placed on the sum over X u . We have 

S(k) = -i- / pe ik »(0\T [Y nd *Yt\^)T[Y na c Yl\m (59) 



using the notation of 

The renormalization of the dijet shape function can be determined from Z 2 . The BPS 
field redefinition |3] is given for all collinear directions n by 

£n * Y^ £ n , 



A — ». yt A' V 



where the path-ordered Wilson line of usoft gluons in the n direction is 



YJz) 



Pexp 



i 9 



dsn ■ A u (ns + z) 



(60) 
(61) 

(62) 



This field redefinition removes all couplings in the Lagrangian of the redefined usoft and 
collinear fields (£' n ,£' n , A' n ). This allows the matrix element to be factorized 

J2 |(J„^X u |xn 9 r^ n , 5 (0)e CT |0)| 2 = J2 \(JnJn\x' n J a x' ncl (0)^m 2 

final states J n Jn 

xJ2\(X u \T[Y n Yt}(0)\0)\*. (63) 



The PSRG is inspired b y th e velocity renormalization group, although some important differences in im 



36 



37 



38 



39( on the equivalence 



plcmcntation exist. See [35| for the VRG in NRQCD. The arguments of 
of one stage and two stage indicate that two stage running could be used to obtain the same results. 
9 However, \X U ) is still not a color singlet. 
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Z°M = 1 + ^%^(4 + f - I In ph.]) + 0{a% (64) 



The renormalization of the initial matrix element must be reproduced by the product of the 
renormalized matrix elements after the field redefinition. Thus Z 2 = Z^Z^^ where Z| is the 
renormalization for the matrix element of only collinear degrees of freedom, while ZJ gives 
the renormalization for the usoft matrix element. These renormalization factors are 

£)Cf ( 2 3 1 v\f 
2tt \ 2 2e e 1 /4 

and 

Z 2 . W = 1 + + I In [1^1) + 0(02). (65) 

Note that the convolution of the variable k between the jet and usoft matrix element is 
neglected as we are considering the total decay rate. The anomalous dimensions of the 
factorized matrix elements are given by 

1 C M = -^(- In [#] + h + 0(a% (66) 

TJ W = -f?£^(l„[!l^i]) +0 K). (67) 

The phase space scales that we will run the usoft and collinear factorized matrix elements 
to are related by = B 5 ~ l /x^/(M|). We relate both fi s and /x c to a phase space parameter 
that runs from 1 to 5 using 

fi 2 s = tfMlB^. (68) 
We can express the running of C2 in the PSRG as 

^C 2 (0) = ( 72 c (0) + (2 + t log(B)) 72 S (0))C 2 (0). (69) 
We take the invariants to be Pg,P? = b 2 M\ and n ■ p q n ■ Pq = f(B) 5 4 M\ giving 



, d a s (M z( p)C F (3 



S 



C^) (70) 



.^j^^^^))^ (71) 

The particular functional dependence on B in the logarithm is fixed by matching onto the 
full theory. Note that we must refer to QCD with the SW jet definition at the collinear scale 
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as we matched the amplitude in QCD onto our SCET operator basis. As this matching was 
performed before the phase space integrals are performed we must refer to the full theory to 
establish the jet definition at the collinear scale. As in SW, we expand in the dependence on 
B and 5 and neglect the dependence on these parameters that does not contribute to large 
logarithms. 

This RGE is solved using standard methods to give 



crM-wm {%ffl*r***"> (72) 

The contributions from the running of the three jet operators contribute at 0(a 2 s ). We 
choose to retain two stage running for the three jet contribution to the Sudakov resumed 
dijet decay rate. 



VII. TWO JET DECAY RATE 



Using C2 (fj/^,fig) we obtain the total decay rate 

* 2 

v(z - j J)(^) = / j£ K0£>IV,) { c 2 W (4^s) + a '^ Cp [\ - 7 -f)){73) 



where we have indicated the dependence on the perturbative expansion of | ) | 2 - The 
contributions to this expression at 0(a s ) come from the running of O2 in the PSRG and the 
mixing of O3 and O2 beginning at the collinear scale. This is sufficient to resum the double 
Sudakov logarithm and the class of sub leading logs we are interested in. The running 
of O3 contributes logarithms that first contribute at 0{a 2 s ). To extend the resummation 
consistently to sub-leading logarithms, the running of O3 must be reexpressed in the PSRG. 
This is beyond the scope of this paper. 

The renormalized perturbative expansion o f K Oq)! 2 is IR finite, as the SW jet definition 



was constructed to satisfy the KLN theorem 



32 



33J. In our effective theory construction, 



the KLN theorem for the SW states ensures that the dijet partial cross section is free of 
explicit IR singularities, as the IR of the full theory is reproduced in the effective theory. 

The expansion of the operator gives the collinear gluon emission contained within the 
dijet cones. However, the fixing of the SW definition within the effective theory removes 
gluons that are excluded by the SW jet definition from the jet cones. This results in constant 
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terms due to the degrees of freedom removed by the jet definition that have to be accounted 
for by matching in the effective theory. We introduce the phase space Wilson coefficient 
S 2 {5 2 Mz) to account for this further matching onto the jet definition in the full theory. As 
a result, we have 

Cp_ i 



CfV„ tf> = C 2 (M Z) 8tf M.) (*gg>) "° (*) . (74) 

We determine S 2 (5 2 Mz) by calculating the SW dijet decay rate in full QCD, and in SCET 
with the PSRG. We match onto the SW jet definition in the full theory to determine the 
perturbative expansion of S 2 (S 2 Mz)- 

Computing the SCET graphs in dimensional regularization we find the following 0(a s ) 
corrections to the matrix element 10 

/ E \{Ol)\ 2 {^) = a {e)(l + 2^^- f 1 dx x [' dx 2 S{x u x 2 , $) A(x u x 2 ))(75) 

J states \ ^IftJ J J(l-xi) j 



with 

5(x 1 ,z 2 ,A*i) 



2 7T r(i - e) (i - Xl y (i - x 2 y ( Xl + x 2 - ly 



A( \ l + x\ e{l-x 2 ) 1 1 ~ e 

{XUX2> (l-x 1 )(l-x 2 ) (l-xx) (l-x 1 )(l-x 2 ) 2{1- Xl )' 1 ' 

Integrating these expressions we find 

/El^)l^)^„ W ( 1 + ^(f-^)). (77) 



states 



To check this result we re-expand the Sudakov factors out into fixed order perturbation 
theory in terms of a s (Mz), we find 

V{Z^JJ){M Z ) = N c M ^ 2 v + 9\) ( 1 + ^M z )C F [ _ 3lnS _ 4HnB)s)l \ 

12 7T \ 7T / 



N c M z (g 2 v + g 2 A ) ( a s (M z )C F 
12vr { ix 



57 7tt 2 

T ~ ~Y 



2a s (M z )s 1 )+0(a 2 s ). 



We match onto the SW dijet decay rate to obtain f(B)S = 2/3. The phase space Wilson 
coefficient S 2 (5 2 Mz) has the perturbative expansion 

S 2 (M Z ) = 1 + a s (M z ) Sl + 0(a 2 s ), (79) 



10 See [29( for details on these effective theory calculations and the requirement of zero bin subtractions. 
Note that as we are integrating over the full phase space. To avoid double counting we find the zero bin 



subtracted contribution to be Rp + Rb — Ra — l/2i?c — 1/2Re in the notation of [29 1 . 
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where 

S ' = c ^-S)- (80) 

This is the main result of our paper. The Sudakov resumed dijet decay rate in SCET is 
given by 

r ( ^ W ) = .„(i + ^(f-^)) (8i) 

* 2 

iSWf 2 2^ , oi s (n 2 s )C 3 C F ( 9 77T 2X 



2tt V 2 6 



Comparing to the literature, we note that a leading log resummation was determined 
in in full QCD for the SW jet definition. The exponentiation of the double Sudakov 
logarithm agrees in both formalisms. Comparing the results, we consider the resummation 
in SCET with the PSRG to be conceptually clearer and easier to extend to non leading 
logarithms. 



VIII. CONCLUSIONS 

We have determined the RGE improved decay rate T(Z — > J J) using SCET. Recalling 
the SW jet definition, we can now state the corresponding definition of the dijet final states 
in our effective theory. 

The dijet decay rate comes from the direct matching of O2 onto QCD and from the mixing 
of the three jet operators at the scale defining the angular cut fij. These contributions 
correspond to the energetic quarks of SW1 and the gluons of collinear energy emitted inside 
the jets, SW3. 

The nonperturbative corrections given by matrix elements of the ultrasoft final state 
gluons are restricted by the scale /a*. These final state gluons correspond to the emission 
of gluons unrestricted in direction with low energy. The ultrasoft matrix element with this 
restriction corresponds to the state SW2. 

Treating the decay rate in SCET allowed the dijet decay amplitude to be defined in terms 
of SCET jet fields. The EFT formalism allowed this definition to take place without large 
logarithms. The fraction of three parton events contributing to the dijet rate in the SW jet 
definition, corresponds to the three jet operators mixing into the dijet operators at the scale 
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fj/g, and the perturbative expansion of 0£ at the scale fig. The general version of this result 
is that n + i jet operators mix into n jet cross sections and decay rates at order a l s (fi 2 6 ). 

To obtain the resummed dijet decay rate, we introduced a one stage running formalism 
for the dijet operator. This allowed us to run to the phase space cut scales defining the dijet 
corresponding to the SW jet definition. We also introduced phase space Wilson coefficients. 
These are required to fix the SW jet definition onto the jet operators in SCET. 

This approach establishes a jet definition in SCET and resums the large phase space 
logarithms that come about as the result of the jet definition. This allows a program of 
systematically improving the perturbative behavior of jet observables to be carried out in 
SCET. 
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□ 

As this paper was approaching completion the work [40J] , which is a more detailed account- 
ing of the work reported on in 30] appeared. In [40} the calculation of the renormalization 



of three jet operators equivalent to the three jet operators examined in this paper, Eq. (1401) . 
was reported. We find that 4(j], which was not calculated in the background field method, 
is inconsistent in its inclusion of the wavefunction renormalization of the gluon while not 
including the QCD coupling constant renormalization. Once 
inconsistency, their results agree with ours. 



30 



is corrected for this 



APPENDIX A: FEYNMAN RULES OF THE 3-JET OPERATORS 

For completeness, we state the one and two gluon Feynman rules of the operators 
03,04,05 to 0(A°Mf) in this Appendix, we utilize these rules to renormalize the oper- 
ators in the Section IIV Al The following Feynman rules are for emitted collinear gluons. 

The one gluon Feynman rules where the gluon is emitted in the direction n 3 requires 
the collinear gluon field strength to emit a gluon. Here P g , P q , Pg are outgoing collinear 



23 



momentum: 

(UUK\O°\0) = U2g s T a 
{UUA» z \Ol\Q) = Ug s T a I 



A 1 



A 1 



(m ■ n 3 )(n 3 ■ P g ) (n 2 • n 3 ) (n 3 ■ P g ) 
1 1 \ 



(ni • n 3 )(ni • P 9 ) (n 2 • n 3 )(n 2 • P ? ) 
1 1 



L Sn2 ) 



(rai • n 3 )(n 1 ■ P q ) (n 2 ■ n 3 )(n 2 ■ Pq) 
The two gluon Feynman rules with P g l and Pg the two outgoing collinear gluon momenta, 



■ 



associated with Ati and Ai are 



2g 2 s (T a T b ^ 



e 



1 fa • Pi) (ns ■ Pg 2 ) 
2gl(T b T a ) 



n 3 n\ — n± n v z n 3 n v 2 — n v z 
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n 2 ■ n 3 
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(n 3 • P, 1 ) (n 3 • P|) 

4^ 2 r fcc T c K<-<^) 
(^•(Pi + Pi)) 2 ^!-^)^-^) 
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L Sn 



r CT f 
r CT £„ 2 

(Al) 



ns ■ PI 
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(n 2 • n 3 )(n 2 ■ P 5 ) 



n 2 



- e 



- g 2 s if abc T c nl (nZTK-n 3 -Tg°») (n° 3 T>* - n 3 ■ T g°> 



n 3 ■ P 2 



(ri! ■ n 3 )(ni ■ P q ) (n 2 ■ n 3 )(n 2 • Pq) 



2g2f"bc Tc / F M gov _ T u g a» g" 7 u - T u g a ^ \ 

+ 1 (n 3 • (Pi + P 2 ) ^(m-ns) (rn-P,) " (n 2 • n 3 ) (n 2 -P,)J 6 
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(finU^^l^lO) = Utfir^T. (ie^M^sT,) - ) (A3) 

\™3 ' -Tg ™3 * } 
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